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.  Application  of  Solidification  Theory  to 
Rapid  Solidification  Processing 
1.  Technical  Report  Summary 

This  semi-annual  technical  report  for  ARPA  Order  3751  covers  the 
period  October  1,  1981  to  March  31,  1982.  Significant  accomplishments 
during  that  period  were  achieved  In  the  following  areas: 


Extended  Solid  Solubilities  --  for  Production  of  Alloys  with  New 
Compositions  and  Phases 

o  At  high  growth  rates  In  Ag-Cu  alloys,  mlcrosegregatlon-free 
single  phase  structures  were  obtained  for  all  compositions 
Investigated,  thus  confirming  that  extended  solid  solubility 
could  be  obtained  completely  across  the  phase  diagram  In  this 
allqy  system. 

o  It  was  found  that  the  solidification  velocity  required  to  produced 
this  single  phase  structure  was  lowest  for  dilute  alloys  and 
alloys  near  the  eutectic  composition,  consistent  with  theoretical 
predictions. 

o  At  Intermediate  growth  rates,  for  Ag-Cu  alloys  with  9X  Cu  or  greater, 
structures  consisting  of  alternate  bands  of  segregated  and 
segregation-free  material  were  obtained. 

Interface  Stability  —  for  Production  of  More  Homogeneous  Alloys 

o  Influence  of  trace  Impurities  on  Interface  stability  during 
solidification  of  dilute  A1  alloys  was  analyzed  and  strong 
destabilizing  effect  for  some  elements  was  found. 

o  Structural  variation  with  composition  which  cannot  be  explained 
In  terms  of  conventional  cellular  patterns  was  observed  in 
In  Al-Mn  alloys. 
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o  Electron  beam  surface  melting  procedures  Mere  developed  which 
allow  the  formation  of  continuous  deep  rapidly  solidified 
layers. 

Interface  Kinetics  --  for  Production  of  Alloys  with  Finer  Segregation 
o  A  theory  was  developed  to  predict  the  effect  of  Interface 

attachment  kinetics  on  the  growth  rate  of  dendrites  In  supercooled 
metals. 

o  A  kinetic  coefficient  equal  to  200  cm/(sK)  was  estimated  for 

pure  nickel  from  the  sound  speed;  this  value  provides  a  good  fit 
to  dendritic  growth  data  at  high  velocities. 

Prevention  of  Cracking  Caused  by  Residual  Stress  --  for  Production  of 
Alloys  with  Improved  Surface  Structures 

o  Response  of  three  different  types  of  metals  to  the  residual 
stresses  developed  during  electron  beam  surface  melting  was 
characterized. 

Task  Objective 

The  objective  of  this  work  Is  to  develop  guidelines  based  on  kinetic 
!  and  thermodynamic  solidification  theory  for  prediction  and  control  of 

rapid  solidification  processes.  In  particular,  segregation  effects  and 
rules  governing  the  formation  of  equilibrium  and  non-equilibrium  phases. 
Including  metallic  glasses,  will  be  Investigated.  Areas  where  significant 
liqprovements  In  allqy  properties  can  be  produced  by  rapid  solidification 
will  be  emphasized. 

Technical  Problem  and  General  Methodology 

Rapid  solidification  techniques  make  It  possible  to  produce  new 
types  of  materials  having  significantly  better  properties  than  conventionally 
processed  materials.  However,  Improved  predictive  techniques  and  control 
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of  rapid  solidification  processes  are  needed.  The  current  studies  are 
focussed  on  the  science  underlying  areas  where  lap roved  materials  can  be 
obtained  In  order  to  provide  such  prediction  and  control.  Ibis  work  Is 
both  theoretical  and  expert mental. 

Three  major  ways  In  which  rapid  solidification  technology  provides 
Improved  materials  are: 

A.  Production  of  alloys  with  new  compositions  and  phases 

B.  Production  of  more  homogeneous  alloys. 

C.  Production  of  alleys  with  improved  microstructures,  Including  finer 
segregation,  and  Improved  surface  structures 

The  general  method  followed  In  this  work  has  been  to  identify  critical 
questions  In  these  three  major  rapid  solidification  application  areas 
where  solidification  theory,  when  properly  developed  and  checked  by 
experiment,  can  provide  Improved  understanding  of  important  rapid 
solidification  processes.  This  understanding  then  is  pursued  to  provide 
guidelines  that  can  be  used  by  alloy  producers  to  obtain  new  Improved 
materials  and  to  select  optimum  alloy  compositions  and  processing  con¬ 
ditions  for  rapid  solidification  applications. 

Accomplishments  to  date  on  this  contract  and  work  In  progress  in 
each  of  these  areas  will  now  be  described  in  more  detail. 

Suwnary  of  Technical  Results,  important  findings.  Plans  and  implications 
for  future  Work 

In  application  A,  particular  success  has  been  achieved  during  the  first 


three  years  of  this  contract  In  (1)  evaluating  conditions  which  control 
critical  solid! ficatfbn  velocities  for  production  of  metallic  glass  alloys. 


and  (2)  determining  conditions  that  control  achievement  of  extended  solid 
solubilities  in  rapidly  solidified  alloys.  The  work  on  critical  solidifi¬ 
cation  velocities  for  metallic  glass  formation  is  nearing  completion. 

Emphasis  now  is  on  investigations  of  extended  solid  solubilities.  Important 
findings  and  planned  future  work  on  these  metallic  glass  and  extended 
solid  solubility  types  are  discussed,  respectively,  in  the  following  two 
paragraphs. 

(1)  By  rapid  quenching  of  Pd-Cu-Si  alloys  of  varying  copper  composition, 
NBS  investigators  were  able  to  determine  maximum  solidification  velocities 
for  dendrite  and  eutectic  crystallization  and,  consequently,  the  conditions 
for  glass  formation.  A  sharp  transition  from  crystalline  to  amorphous 
solidification  was  consistently  found  when  the  solidification  velocity 
was  increased  above  about  2-20  mm/s  in  these  alloys,  depending  on 
composition.  These  results  then  were  used  to  relate  glass-forming 
tendencies  to  phase  diagram  features  as  well  as  to  kinetic  models  of 
crystallization  based  on  the  diffusional  sorting  of  solute  in  the  liquid. 

At  some  compositions  and  velocities,  independent  nucleation  of  crystallites 
occurred  even  in  the  presence  of  a  substrate.  Nevertheless,  the  experimental 
results  generally  confirmed  the  importance  of  solute  redistribution  on  the 
maximum  crystal  growth  rate  in  these  alloys.  Crystal  growth  rates  rather 
than  nucleation  rates  were  the  main  factor  in  determining  whether  amorphous 
or  crystalline  solidification  was  obtained.  In  carrying  out  this  work,  a 
directional  solidification  apparatus  was  built  which  allowed  direct 
knowledge  of  fast  solidification  velocities.  These  velocities  usually 
have  not  really  been  known  in  previous  rapid  solidification  investigations. 
(2)  Principles  which  allow  prediction  of  maximum  thermodynamic  limits  on 
extended  solid  solubility  in  rapidly  solidified  alloys  and  phase  sequences 
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for  metastable  alloys  have  been  formulated.  However,  because  of  kinetic 
limitations,  it  may  in  some  cases  not  be  possible  to  obtain  the  amounts 
of  extended  solid  solubility  that  would  be  thermodynamically  allowed. 
Experiments  are  being  performed  to  test  the  predicted  kinetic  and 
thermodynamic  limits  on  extended  solid  solubility.  Here  again,  it  will 
be  important  to  establish  the  rapid  solidification  velocities  that  are 
sufficient  to  prevent  the  sorting  of  constituents  into  two-phase  rather 
than  one-phase  structures.  Electron  beam  surface  melting  and  refreezing 
of  Ag-Cu  alloys  is  being  used  to  investigate  these  effects  more  thoroughly. 

In  application  B,  it  was  predicted  that,  when  solidification  occurs 
rapidly  enough,  the  stable  shape  of  the  solidification  interface  is  planar 
rather  than  cellular  or  dendritic.  Then,  interdendritic  alloy  segregation 
will  be  avoided  and  a  more  nearly  homogeneous  alloy  obtained.  Planar 
interfaces  are  also  stable  at  very  slow  velocities.  However,  at  intermediate 
velocities,  nonplanar  (cellular  or  dendritic)  interfaces  are  produced. 
Calculations  indicate  that  the  high  velocity  planar  regime  will  occur  at 
practical  rapid  solidification  rates  in  a  number  of  alloys.  Experiments 
are  being  performed  on  aluminum  alloys  to  test  this  prediction.  An 
electron  beam  melting  system  has  been  modified  to  produce  a  measurable 
energy  distribution  in  the  beam  during  surface  melting.  When  this  beam 
is  pulsed  or  scanned  across  the  surface,  solidification  velocities  produced 
during  the  refreezing  process  can  be  calculated  to  allow  comparison  of 
experimental  observations  with  theory.  Conditions  under  which  stable  inter¬ 
faces  can  be  maintained  are  being  investigated  as  a  function  of  solidification 
velocity  and  alloy  composition.  While  elements  having  high  solute  dis¬ 
tribution  coefficients,  k,  such  as  manganese  and  silver,  should  not  be 
destabilizing  even  at  rather  high  concentrations  (fractions  of  one  percent) 
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at  solidification  velocities  of  1  to  10  cm/s,  results  to  date  Indicate  that 
other  solutes,  such  as  nickel,  tin,  or  Iron,  having  much  lower  distribution 
coefficients  can  cause  instability  at  these  same  velocities  when  present 
In  much  lower  concentrations.  Investigations  of  the  microstructures  formed 
in  these  alloys,  their  causes,  possible  effects  of  non-equilibrium  at  the 
Interface,  and  implications  for  rapid  solidification  processing  are  continuing. 

In  application  C,  three  effects  of  rapid  solidification  on  alloy 
microstructure  are  being  investigated;  (1)  microsegregation  in  rapidly 
solidified  alloys,  (2)  formation  of  amorphous  layers  on  alloy  surfaces, 
and  (3)  stresses  leading  to  cracking  of  rapidly  solidified  surfaces  are 
being  investigated. 

(1)  A  large  number  of  alloys  can  be  expected  to  exhibit  cellular 
growth  and  consequent  fine  scale  intercellular  segregation  even  at  very 
high  solidification  velocities.  To  investigate  this  segregation  a  quanti¬ 
tative  theory  of  microsegregation  during  rapid  solidification  has  been 
formulated  which  shows  the  influence  of  cellular  spacing,  solidification 
velocity,  and  the  amplitude  of  the  instability  at  the  interface  on  the 
magnitude  and  distribution  of  the  intercellular  segregation.  This 
fine-scale  segregation  is  important  because  of  the  precipitation  processes 
that  may  be  produced  and  because  of  the  well-known  correlation  between 
segregate  spacing  and  mechanical  properties. 

(2)  Methods  have  been  devised  for  producing,  on  Pd-Cu-Si  substrates, 
by  electron  beam  surface  melting  and  refreezing,  continuous  amorphous 
layers  up  to  200  pm  thick  without  optically  visible  crystals  being 
formed.  The  formation  of  glassy  surface  layers,  particularly  if  extended 
to  other  alloys,  could  produce  special  desirable  surface  properties  on 
otherwise  crystalline  materials.  The  overlapping  strip  method  used  to 
produce  these  deep  glassy  layers  is  also  being  applied  to  studies  of 
surface  layers  in  crystalline  materials. 
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(3)  Several  deleterious  effects  can  result  from  the  residual 
stresses  generated  during  surface  melting,  such  as  cracking  in  the  melt 
zone  or  substrate  and  distortion  of  the  substrate.  Cracks  produced  after 
electron  beam  surface  melting  and  refreezing  at  various  sweep  velocities 
and  power  inputs  are  being  investigated  to  provide  possible  guidlines  for 
avoidance  of  this  persistent  problem  which  occurs  in  many  instances  where 
rapid  solidification  is  used  for  surface  modification  purposes.  Materials 
with  less  extreme  susceptibility  to  plastic  deformation,  hot  tearing  and 
brittle  fracture  will  be  used  to  define  more  clearly  the  variation  of 
defect  structure  with  processing  parameters. 


2.  Report  of  Technical  Progress  and  Results 
Extended  SoTid  Solubility 

The  goal  of  this  research  is  to  determine  experimentally  the  solidification 
velocity  required  to  produce  microsegregation-free  single  phase  alloys  in 
the  Ag-Cu  alloy  system.  For  alloys  exceeding  the  equilibrium  solubility, 
such  a  microstructure  indicates  extended  solubility.  The  Ag-Cu  system 
has  been  chosen  because  the  phase  diagram  has  a  simple  eutectic  between 
two  fee  phases  and  the  thermodynamic  functions  are  sufficiently  well- 
known  to  permit  a  determination  of  the  T0  curve,  which  is  the  locus  of 
points  on  the  phase  diagram  where  the  liquid  and  solid  free  energies  have 
the  same  value.  In  addition,  the  metastable  solidus  curves  of  the  Ag- 
and  Cu-rich  phases  appear  to  be  retrograde  such  that  solubility  extension 
beyond  these  retrogrades  must  be  accomplished  by  solute  trapping. 

Two  solidification  mechanisms  can  produce  segregation-free  crystalline 
solids:  planar  growth  and  partitionless  solidification.  The  solidification 
conditions  required  for  planar  growth  with  equilibrium  partitioning  at 
the  interface  are  well  known.  For  growth  at  high  velocities,  capillarity 
provides  a  stabilization  of  the  planar  interface.  This  regime  is  referred 
to  as  absolute  stability  and  applies  only  when  the  net  heat  flow  is  toward 
the  solid.  Partitionless  solidification  arises  from  the  kinetics  of 
interface  motion.  Solute  is  trapped  by  the  moving  interface  causing  the 
partition  coefficient  to  deviate  from  the  equilibrium  value.  Solute 
trapping  occurs  at  high  interface  velocities  when  the  interface  temperature 
Is  below  a  temperature  determined  from  thermodynamic  considerations.  In 
the  limit  of  partitionless  solidification,  when  the  partition  coefficient 
is  unity,  this  temperature  is  the  T0  temperature  for  the  composition  at 
the  Interface.  A  major  goal  of  this  research  Is  to  determine  the  relative 
contribution  of  these  two  mechanisms. 
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Alloys  were  prepared  for  electron-beam  surface  melting  by  induction 
melting  of  99.99%  pure  components  in  a  graphite  crucible  under  vacuum  and 
chill  casting  in  a  carbon-coated  copper  slab  mold.  Ingots  were  cold- 
rolled  to  3  mm  thick  sheets.  For  the  surface  melting  experiments,  the 
substrate  must  be  prepared  with  either  no  microsegregation  or  microsegregation 
fine  enough  to  ensure  complete  mixing  during  surface  melting.  The 
following  procedure  was  found  to  be  adequate  for  solidification  at 
velocities  up  to  400  cm/s.  Alloys  containing  1  and  5  wt%  Cu  were 
homogenized  for  20  hrs.  at  750  °C.  Alloys  containing  9,  15  and  23  wt%  Cu 
were  surface  melted  using  overlapping  electron-beam  passes  to  a  depth 
of  ~  0.5  mm  at  a  scan  speed  of  3  cm/s.  This  produces  a  dendrite  arm 
spacing  in  the  range  of  1  pm.  These  slabs  were  then  rolled  slightly  to 
eliminate  surface  topography  produced  by  the  e-beam  melting.  Alloys 
containing  28.1  wt%  Cu  (the  eutectic  composition)  required  no  pretreatment 
due  the  fineness  of  the  eutectic  structure. 

Surface  melting  experiments  were  performed  by  electronically  scanning 
a  25  keV  electron  beam  focussed  to  approximately  1  mm  with  a  current  of 
50  mA  across  the  surface  of  the  prepared  substrates.  Scan  ve.ocities 
between  1.5  cm/s  and  400  cm/s  were  possible.  At  the  slow  scan  speeds 
welds  were  typically  0.6  mm  wide  and  0.2  mm  deep  while  at  400  cm/s  the 
welds  were  typically  0.2  mm  wide  and  0.02  mm  deep.  Samples  for 
metallography  weVe  obtained  as  sections  transverse  to  the  scan  direction 
and  occasionally  the  other  two  perpendicular  sections.  TEM  samples  were 
obtained  by  punching  half  disks,  3  mm  in  diameter,  from  0.1  mm  thick 
transverse  slices  through  the  weld  zone.  Thin  foils  were  obtained  by  ion 
milling. 
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Significant  effects  arise  from  the  change  in  weld  puddle  shape  as 
the  scan  rate  is  Increased.  The  local  solidification  velocity  is  equal 
to  the  product  of  the  scan  velocity  with  the  cosine  of  the  angle  between 
the  scan  direction  and  the  local  growth  direction.  Because  of  this 
factor,  growth  velocities  at  the  top  center  of  the  weld  puddle  for  high 
scan  speeds  may  differ  from  the  scan  velocity  by  a  factor  of  two.  Four 
general  types  of  microstructure  are  observed  depending  on  solidification 
velocity  and  composition,  a)  At  high  growth  rates,  microsegregation-free 
single  phase  structures  are  obtained  for  all  compositions.  The  velocity 
required  to  produce  this  structure  is  lowest  for  dilute  alloys  and  alloys 
near  the  eutectic  composition.  At  low  growth  rates,  b)  conventional 
dendritic  (cellular)  or  c)  eutectic  structures  are  obtained,  d)  At 
intermediate  growth  rates  for  alloys  with  9  wt%  Cu  or  greater,  a  structure 
consisting  of  alternate  bands  of  segregated  and  microsegregation-free  material 
is  obtained.  The  bands  form  parallel  to  the  local  interface.  The  growth 
rates  required  to  cause  the  various  changes  in  microstructure  for  the 
different  compositions  will  be  described  in  the  future. 
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Experiments  on  Interface  Stability  and  Segregation  with  Electron  Beam  System 

Electron  beam  surface  melting  can  provide  a  very  well -control led 
method  for  producing  rapidly  solidified  regions  on  the  surface  of  bulk 
metal  substrates.  The  electron  beam  facilities  at  NBS,  and  the  techniques 
for  operating  them,  have  been  continuously  upgraded  so  that  it  has  become 
possible  to  better  characterize  the  melting  conditions  and  to  make  deeper, 
smoother,  and  more  uniform  layers  of  rapidly  solidified  material  for 
metallurgical  evaluation. 

One  important  objective  of  rapid  surface  solidification  is  to  produce 
highly  homogeneous  surface  layers.  We  report  here  work  relating  to 
achieving  this  objective  by  two  different  methods; 

(a)  by  crystalline  solidification  at  high  velocities  and 

(b)  by  the  formation  of  continuous  layers  of  metallic  glass 

Analysis  of  the  influence  of  several  different  solutes  on  the 

morphological  stability  of  aluminum  during  rapid  solidification  has 
indicated  that  while  elements  having  high  solute  distribution  coefficients, 
k,  such  as  manganese  and  silver,  can  be  present  in  rather  high  concentrations 
(fractions  of  one  percent)  without  causing  interfacial  instability  at 
solidification  velocities  of  1  to  10  cm/s,  other  solutes  such  as  nickel, 
tin  or  iron  having  much  lower  distribution  coefficients  can  cause 
Instability  at  these  same  velocities  when  present  in  much  lower 
concentrations,  on  the  order  of  10"^  percent.  Assuming  equilibrium 
partitioning  of  solute  at  the  solid-liquid  interface,  the  relative 
Influence  of  the  low  k  elements  is  greater  at  higher  velocity  (the  absolute 
stability  regime).  If,  however,  the  solidification  velocity  is  sufficiently 
high  to  give  non-equilibrium  partitioning  effects  (solute  trapping),  then 
the  relative  influence  of  low  k  elements  may  be  greatly  diminished.  In 
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the  velocity  region  where  the  distribution  coefficient  varies  strongly 
with  velocity,  the  stability  analysis  becomes  much  more  complex. 
Unfortunately  very  little  information  Is  available  concerning  the  velocity 
dependence  of  k. 

As  a  result  of  these  considerations  it  is  probable  that  minute 
amounts  of  impurities  such  as  iron  account  for  our  observation  that 
aluminum  containing  Mn  or  Ag  is  less  stable  than  predicted  at  solidification 
velocities  of  a  few  centimeters  per  second  but  can  nonetheless  be  made 
stable  at  high  concentrations  by  solidifying  at  velocities  which  may  be 
sufficient  to  produce  solute  trapping  (as  produced  by  melt  spinning). 

Metal lographic  analysis  of  electron  beam  melt  spots  on  Al-Mn  alloys 
showed  a  structural  variation  with  composition  which  is  difficult  to 
explain  in  terms  of  a  conventional  cellular  structure  and  is  suggestive 
of  velocity-dependent  distribution  coefficient  effects.  Conditions 
which  in  a  1%  Mn  alloy  produce  a  cellular  microstructure  which  is  very 
strongly  lamellar  in  a  longitudinal  section  produced  in  an  0.1%  Mn  alloy 
a  structure  which  in  the  same  longitudinal  section  appears  much  more  like 
a  three-dimensional  network  of  cells.  A  complete  characterization  of 
this  structure,  however,  will  require  simultaneous  determination  of  its 
appearance  in  orthogonal  sections.  Comparison  with  the  theory  of  interface 
stability  including  the  effects  of  a  velocity-dependent  distribution 
coefficient  is  rendered  difficult  by  the  absence  of  any  information  on 
the  form  of  k(v)  and  the  possible  influence  of  minute  concentrations  of 
Impurities.  ' 

The  most  extreme  case  of  a  homogeneous  metallic  solid  is  the  metallic 
glass,  and  we  have  found  that  electron  beam  surface  melting  can  be  used 
to  make  relatively  thick  extensive  layers  of  amorphous  materials.  This 
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study  used  palladium-copper-silicon  alloys  in  which  it  is  known  that 
coupled  growth  of  the  crystalline  phases  cannot  take  place  at  velocities 
of  greater  than  a  few  millimeters  per  second.  Such  sluggish  growth  kinetics 
suggested  the  possibility  that  if  the  electron  beam  is  scanned  rapidly 
across  the  surface  to  produce  a  single  amorphous  melt  pass,  which  is 
known  to  be  quite  easy  in  this  material,  then  successive  overlapping 
passes  could  be  used  to  generate  additional  amorphous  material  without 
resulting  in  significant  crystallization  in  the  previous  passes.  Thus 
even  though  the  amorphous  material  is  raised  to  a  high  temperature  by 
successive  passages  of  the  electron  beam,  any  crystals  which  might  form 
cannot  grow  more  than  a  few  microns  before  the  local  region  again  cools  to 
ambient  temperature. 

The  strategy  used  to  produce  continuous  amorphous  layers  was  thus  to 
start  with  a  crystalline  slab  6-8  mm  thick,  form  a  line  of  amorphous 
material  by  a  single  sweep  of  the  electron  beam  across  the  surface,  allow 
sufficient  time  for  the  sample  to  cool  thoroughly,  and  then  sweep  the 
electron  beam  across  the  surface  again,  offset  from  but  overlapping  the 
I  first  sweep.  The  process  was  continued  to  generate  a  strip  of  any  desired 

width  of  overlapping  melt  passes.  The  time  allowed  between  successive 
passes  was  30  or  60  seconds.  This  time  was  more  than  enough  to  allow  the 
amorphous  material  to  cool  to  a  temperature  close  to  that  of  the  substrate, 
but  was  kept  long  to  allow  sufficient  transfer  of  heat  from  the  substrate 
to  the  water-cooled  copper  chill  plate  below  it.  The  Individual  melt 
passes  were  approximately  1  mm  wide  and  the  offset  between  passes  was 
typically  0.5  mm. 

Continuous  amorphous  strips  were  formed  by  this  process,  although 
usually  with  extensive  cracking  of  the  substrate  as  discussed  in  another 
part  of  this  report.  Crystals  which  were  present  within  the  melt  zones 
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had  microstructures  Indicative  of  growth  radially  outward  from  nucleatlon 
sites,  thus  Indicating  that  the  sweep  speed  of  the  electron  beam  Is 
significant  not  because  of  any  possibility  of  the  crystals  growing  with 
this  speed  but  because  It  determines  the  time  which  the  local  region  may 
spend  at  elevated  temperature. 

Crystallization  of  previously  formed  metallic  glass  by  subsequent 
melt  passes  was  concentrated  near  the  bottom  of  the  melts  with  the  result 
that  a  continuous  layer  of  glass  remained  near  the  surface  if  the  sweep 
speed  was  sufficiently  high.  Material  transport  by  the  moving  beam  was 
extremely  prominent  in  these  experiments,  with  the  result  that  the  final 
surface  was  greatly  depressed  below  the  original  surface  in  the  region 
where  the  melt  passes  started,  and  elevated  above  the  original  surface  by 
a  corresponding  amount  at  the  opposite  end. 

Using  the  overlapping  scan  method,  continuous  glassy  layers  up  to 
200  pm  thick  without  optically  visible  crystals  were  formed.  This 
required  sweeping  the  electron  beam  across  the  surface  at  velocities  of 
60  cm/sec  or  more.  Further  optimization  of  processing  parameters  could 
result  in  thicker  amorphous  layers,  since  the  overlap  method  has  been 
shown  to  work  and  some  glass  was  formed  at  depths  of  over  500  un. 

When  the  same  process  Is  attempted  with  Cu-Zr,  very  little  glass  is 
formed.  Although  this  is  also  a  relatively  easy  glass  former,  when  the 
melts  are  several  hundred  micrometers  deep  the  cooling  rates  are  not 
sufficient  to  prevent  the  crystals  In  the  substrate  from  growing  up  to 
the  surface  in  most  places. 

In  addition,  attempts  to  produce  a  wide  strip  of  thick  glassy  material 
on  Pd-Cu-SI  by  rapidly  (several  kHz)  oscillating  the  electron  beam  in  the 
direction  perpendicular  to  that  with  which  it  swept  over  the  surface  of 
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the  slab  resulted  in  amorphous  material  only  along  the  sides  of  the  strip. 
Although  no  heat  flow  analysis  of  this  melting  mode  has  been  carried  out, 
ft  Is  clear  that  the  attainable  cooling  rates  are  slow  except  at  the 
edges  of  the  strip  or  for  very  thin  melts.  The  overlapping  strip  method 
as  described  above  Is  thus  the  most  effective  way  to  obtain  deep  layers 
of  rapidly  solidified  material,  and  It  is  being  applied  to  studies  of 
surface  layers  in  crystalline  materials,  where  the  transient  heating  of 
the  earlier  melts  is  not  as  critical  as  In  the  glassy  materials. 

Theory  of  Interface  Kinetic  Effects  and  Interface  Stability 

Significant  departures  fro m  local  equilibrium  at  the  crystal-melt  Interface 
can  exist  for  solidification  at  rapid  rates.  Perhaps,  the  highest  measured 
crystal  growth  velocities  (50  m/s)  are  those  of  Ni  dendrites  growing  into 
pure  undercooled  melts  as  reported  by  Walker  and  Colllgan  and  Bayles. 

Recent  theoretical  (Langer  and  Mul ler-Krunbhaar)  and  experimental  (Glicksnan 
et  al.)  advances  In  our  understanding  of  dendritic  growth  allow  accurate 
calculation  of  the  heat  transport  limited  rate  of  dendritic  growth  as  a 
function  of  melt  undercooling.  Inclusion  of  interface  kinetics  into  the 
heat  transport  model  of  dendritic  growth  has  been  performed  in  collaboration 
with  Professor  David  Turnbull  to  extend  the  model  to  high  solidification 
velocities.  We  have  used  previous  results  on  the  effect  of  interface 
kinetics  on  the  morphological  stability  of  a  solidifying  sphere  to 
calculate  the  effect  of  interface  kinetics  on  the  dendritic  tip  radius 
and  hence  the  growth  rates.  An  Interface  kinetic  coefficient  based  on 
the  estimated  velocity  of  sound  In  molten  nickel  provides  a  good  fit  to 
the  dendritic  growth  data  at  high  velocities  as  a  function  of  melt 
undercooling. 
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The  details  of  this  research  are  given  In  an  appendix  In  a  Manuscript, 
"Relative  Reles  of  Heat  Transport  and  Interface  Rearrangement  Rates  In 
the  Rapid  Growth  of  Crystals  In  Undercooled  Melts,"  by  S.  R.  Corlell  and 
0.  Turnbull  to  be  published  In  Acta  Metal lurglca. 

Prevention  of  Cracking  Produced  by  Electron  Beam  Surface  Melting  and  Refreezlnq 
Surface  layers  of  rapidly  solidified  material  can  be  produced  on 
bulk  metal  substrates  by  pulsed  or  scanned  energy  beams  (lasers  or  electron 
beams).  Such  layers  may  have  Improved  hardness,  wear  or  corrosion 
resistance,  or  other  desirable  properties  but  In  many  cases  the  potential 
benefits  are  offset  by  the  harmful  effects  of  the  residual  stresses 
developed  during  the  surface  melting  process.  The  origin  and  the  effects 
of  these  stresses  are  essentially  the  same  as  those  developed  during 
conventional  welding  processes,  but  the  magnitude  of  some  of  the  controlling 
parameters  such  as  temperature  gradients  may  be  significantly  different. 

|  Several  deleterious  effects  can  result  from  the  residual  stresses 

generated  during  surface  melting.  The  most  Important  can  be  listed  as 
follows: 

|  (a)  The  most  severe  type  of  defect  arising  from  surface  melting  Is 

cracking,  which  may  occur  In  the  melt  zone  or  the  substrate 
or  both.  The  prevalence  of  this  type  of  defect  Is  dramatized 
by  the  frequent  publication  of  photomicrographs  of  surface 
melts  containing  conspicuous  cracks,  when  the  authors  would 
presumably  prefer  to  have  obtained  melts  without  cracks. 

(b)  The  residual  stresses  from  surface  melting  can  also  cause  severe 
distortion  of  the  substrate.  This  problem  becomes  dramatically 
clear  when  one  attempts  to  carry  out  surface  melting  of  large 
areas  on  relatively  thin  (less  than  about  5  mm  thick)  substrates. 
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(c)  The  Inhomogeneous  nature  of  the  residual  stresses  can  be  expected 
to  provide  a  altitude  of  favorable  sites  for  stress  corrosion 


cracking  and  other  fores  of  electrochemical  degradation. 

(d)  The  fact  that  surface  melting  almost  always  results  in  a  surface 
which  is  in  a  state  of  tensile  stress  means  that  even  if  cracks 
are  not  present  Immediately  after  solidification,  the  surface 
will  still  be  a  favorable  site  for  crack  initiation. 

The  importance  of  these  different  effects  will  depend  upon  the  alloy 
being  melted,  the  application  for  which  it  is  intended,  and  the  mode  of 
surface  melting  which  is  employed.  Because  little  attention  has  until 
now  been  devoted  to  these  effects,  we  have  started  studies  of  residual 
stress  effects  with  the  ultimate  objective  of  identifying  procedures 
which  may  ameliorate  or  eJImlnate  the  problems. 

The  most  practical  way  to  produce  extensive  rapidly  solidified  surface 
layers  on  a  metal  Is  to  scan  an  energy  beam  over  the  surface  of  a  material 
to  produce  a  continuum  of  overlapping  small  weld  lines.  This  mode  of 
melting  can  be  contrasted  to  the  mode  frequently  used  for  laser  annealing 
of  semiconductors,  In  which  a  pulsed  taser  beam  melts  a  wide  but  extremely 
shallow  layer  on  the  material  surface,  tn  the  latter  case,  the  heat  flow 
Is  one-dimensional  and  residual  stresses  may  be  small  or  absent  but  this 
mode  of  melting  Is  Impractical  for  most  metallurgical  applications  where 
deeper  melts  are  desired. 
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When  a  mil  at  It  spot  Is  scanned  across  the  surf  act  of  a  aetal,  the 
deleterious  offsets  mentioned  above  have  their  origin  in  tht  thermal 
contraction  of  tht  solid  aatarlal  as  It  cools,  and  In  sons  casts  also  In 
tht  change  of  voluat  which  accoapanfts  solidification.  As  tht  scanning 
velocity  of  tht  melt  spot  increases,  the  heated  area  surrounding  the  nett 
zone  changes  In  shape  and  size  and  the  response  of  the  system  to  tht 
stresses  developed  during  cooling  also  changes.  Realistic  analysis  of 
the  stresses  developed  during  this  type  of  process  is  rendered  difficult 
by  the  temperature  dependent  elastlc/plastlc  behavior  and,  especially  In 
the  cast  of  highly  alloyed  materials,  the  possibility  of  large  deviations 
from  egul librium  during  solidification.  Of  particular  importance  in 
surface  melting  of  many  alloys  are  hot  tearing  phenomena  in  trtiich  cracks 
form  by  the  pulling  apart  of  the  material  In  the  solute-rich  regions 
between  dendrites  or  between  grains. 

For  the  case  of  melting  by  a  circular  energy  spot  scanned  across  a  metal 
surface.  It  has  been  shown  (1)  that  the  critical  parameters  controlling 
the  melting  are  qa  and  Ua/2n,  where  q  is  the  absorbed  power  density, 
a  is  the  spot  radius,  U  is  the  velocity  with  which  the  spot  moves,  and  a  is 
the  thermal  dlffuslvlty.  For  a  stationary  spot  on  a  given  material,  a 
critical  value  of  qa  is  needed  to  initiate  surface  melting  and  the  surface 
reaches  the  vaporization  tewperature  at  a  higher  critical  value  of  qa  (2). 

The  latter  condition  bears  a  relation  to  the  development  of  the  deep 
penetration  mode  of  welding,  which  Is  not  generally  desired  for  surface 
melting.  As  a  result  the  range  of  useful  values  of  qa  for  surface  melting 
is  relatively  narrow,  from  approximately  2.3  x  10*  w/m  to  4  x  10*  M/m  for 
aluminum.  When  the  energy  spot  moves,  the  melt  pool  does  not  differ 


greatly  In  ihipt  from  the  ititiaMry  ipot  for  U o/?«  values  1m  Um 
unity.  Although  the  shop*  of  the  atlt  pool  Hot  not  boon  computed  for 
Ue/?«>1,  ft  ft  door  that  tHo  melt  pool  trails  increasingly  for  behind 
tHo  Moving  tpot  and  that  tHo  roor  edge  becomes  increasingly  pointed.  bith 
Increased  velocity  tHo  no  It  pool  olto  diminishes  »n  two  ond  eventually 
dftoppoort  unlott  tHo  qo  voluo  ft  t imultoneoutly  incrootod.  YHoto  cHongot 
tn  tHopo  of  tHo  no U  pool  with  Inc rooting  velocity  lood  to  cHongot  in  the 
rot  I  duo)  ttrott  pottomt. 

Electron  boon  surface  notltng  by  l hit  typo  of  scanned  circular  tpot 
prodocot  dtfforont  rot f duo I  ttrott  rooctfont  In  dffforont  alloys,  ot  ft 
o  spot  tod  on  tHo  botft  of  otoilfc/ploftic  proportiot  ond  solidification 
pHononono.  tn  tono  cotot  tHo  rot I duo 1  ttrott  effects  tory  tfgnff fcontly 
with  tHo  noil  tcon  velocity.  In  our  experiments,  noltt  were  corrfod  out 
by  electronic  deflection  of  o  »»  electron  boon  ot  velocities  from  ?  to 
100  cm/sec.  THo  froo  turfoco  ond /or  potffHod  toctfont  «e  r»  e«*mintd  for 
ovfdonco  of  rot 1 duo  I  ttrott  rolotod  dofoctt. 

tn  puro  otuntnun  crocks  connot  bo  gonorotod  by  oHHor  Hot  toorfng 
during  to) Id) ff cot  ton  or  subsequent  brtttlo  frocturo.  Instead,  thft 
notorfo)  rotpondt  to  tHo  gonorotod  ttrottot  by  plottfc  deformation,  ond 
If  tHo  turfoco  ft  tufffcfontly  clean  on  obundonco  of  ttfp  If net  become 
vftfblo.  THote  tlfp  If net  extend  ocrott  tHo  turfoco  of  the  melt  rone  ond 
extend  without  dftcontfnuftjr  into  the  matrix  on  either  tide.  At  High 
velocities  the  tlfp  lines  ore  more  prominent  ond  ore  visible  ot  o 
larger  distance  to  the  side  of  the  melt  path  than  they  are  at  lower 
velocities.  It  is  not  clear,  however.  If  this  difference  Indicates  a 
greater  amount  of  slip  at  high  velocities,  or  o  concentration  of  the  slip 


toto  i  «w)ltr  amitr  of  —  visible  te«4t.  f**the* 
iii*f  ftoflf  (ryflil  iuMlntM  — y  rtMitt  |li(  owe*  do*. 

hhen  *epe4ted  o—*t  4991*9  —Its  *r»  urritf  owl  o*  pw*0  4t«*«nw*  to 
4  *4pldly  solidified  s*rf4<e  ityfr,  too  froo  «wrf«co  u  co— d 
otto  4  potto**  •#  intersect  1*9  slip  Ii«h  tiid  s(*o*9ly  —lo  to* 
stretches  eo  o  port  lolly  polished  — tottoorophic  I*  so—  coses 

ton#  flip  linos  <00)41  00  too  limtlit  f4<t0*  With  rotfMKt  to  to*  PDOtOOU 
Of  too  treoted  temple  («rfMO,  Oot  —  tty  totf  4—  *41—*  wwll  in 
omplltwde  (Mportt  to  too  0*0f4l  I  twrfoce  roliff  (OwSOd  Of  too  —It  1*9 
process. 

lo  cent rest  to  pwre  4  loot  now,  J03N  4)  lor  oi|0  it*  hip*  cooo or  content 
Show*  0*1 0**1  —  (rocllut  lo  *0*000*0  to  *w*f4<0  —111*9.  A  *o*lo*  of 
oolt*  wo*  00O*  00  101*  otto r  W*l*9  04  *4 two*  >01— OW  )  •  10*  end  *  *  10* 
w/w  4t  — tocltlo*  f*on  ?.S  •  10**  to  I  01/*.  (0**0*00*41*4  to  U4/?W  fee n 
O.OS  to  ?.  In  INt*#  0*00*1— nt*  the  otott*—  >000  *441**  4  «4t  —10 
con* toot  ot  400*0*1— toty  J.S  *  10*4  ft.  HUM*  t>l*  *0*90  of  oo*o— to**, 
t>0  <*0(4109  polio**  »AO—0  tttlto  0*140—0  of  04*14110*  with  04  —1  4 
(too*  —Motion  wlt>  —  toe  tty.  St  Ip  tl— *  — *0  — — *  *oo*  1*  this  —  tori#!. 

Alt  of  l—  C*4Cft*  —It*  — *«  ObtO*— d  40004—9  tO  be  *0t  t«4*«*9  * 
l.o.,  they  occur  rod  In  tolvlewlch  ro^jlo**,  either  Pot— on  9*41**  or 
between  9on9*lto*. 

At  —tocltlo*  of  S  *  1(T*  m/%  0*  to**,  crock*  — ro  — *y  t— 11  end 
—  frequently  Iftterdendrltlc  4*  —11  4*  1*to*9*4—  lor.  they  showed  4 
*tro*9  tendency  to  lie  perpendlcwlor  to  the  loco)  tolid'H—ld  Interfoce, 
the  shopo  of  which  1*  del  ineoted  by  wrfKe  rtoole*. 
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m  a  ttfMt  mImiIj  o#  l<M  m/%.  ititft  omoms  e*i#a*t. 

fNIt  K(»rf  ootl  ly  »•  II*  for#  mi  a  ll*|t*  lirft  (r*(l  ruMitf  |l*»f  !Nt 
(**t*f  I  III*  •#  IN#  *#ll  «#ft,  (r*a«  «r*  *lu  y*wK  r#wt«f 

Mw*N  fr*N  IN#  (<*H  of  IN#  0#l|.  N«(  (My  *f*  wUm  (N#  (M 
«ra<N  «M  *r*  *Mlly  «#f  *ir#<|ly  (NM#(IK  I#  M.  H  iMt  #l*(My  *M 
#|  III  IIUMf  mImIIIH.  IN#  (racNt  «Ui»r  |I«MI  oatlotlaoly  t>«H  «r*t« 

**  IN#  #lM>ly  *M»  «il(l  IN#  electro*  ***#  fQNNt  *#r  (N*  mHic* 
It  iMfNM  «#  I#  I  #/t,  IN#  cram  r«Mi«|  1*  fro*  (M  t*#et  of  IN#  WH 
N#CflNO  «#t  #»##M##|  #*#  IN#  (rMI  |l#|  IN#  C##|#r  ll*»  of  |N#  (#l| 

#Ott  NOCOMOt  r#t«|l«#ly  «M#lt#r  404  4009  «itCOH  i«*0«t.  H  IN#  IlfMU 
t####t.  IN#  lt»#tmy  of  IN#  (r«M  I#  440  f#*»##i><#»#r  i#  (•#  local 
to)l#>)l#ot#  l*l#*ti(#.  #t  i#ti#t#l#i  fey  IN#  torfac#  r«p#l#t,  «#cr#at#t. 

I#  IN#  tolocHy  rg*f#  cl#t#  I#  to*1  mis.  409*9  C*#CN>««#  «u»n 
trloorlly  •*«•#  IN#  li#  of  in#  #ll,  a*  *«f(NM  area  #f  »*l*H#1y 

CfacN«fr#0  t «rf«#  «#l|  torn  t#  #*##**##  ly  O*#rta#0*«#  t#fet#Qei*t  «#t| 
yattat  to*  0009/0  I#  fi*#H  IN#  (p|#r  H«#  cracftt. 

A  totally  <lfr«ffa|  lytf  of  cr#<N <•#  ms  tee*  «•  IN# 
yalta#t«»>coifor  •title#*  alloyt  4000  #«#*  I  «##«•§  «#1t  tract  %  «#r#  «t#a  l# 
KO*K*  «  CMtlMNrt  flatty  torface  layer.  IN#  flatty  layer  4*ct11e, 
oo4  all  cracMftf  otufffl  I#  IN#  tofettrale.  Ir###c#tfy  IN#  feature  t«# 
IN#  tom  of  a  *#ry  lar#t  t«#fl#  crac*,  ly**f  oMfc  IN#  flatty  layer  omo 
at  a  #»#IN  a#  to  two  or  tNre#  l feet  IN#  #c#IN  of  IN#  a# It  rone.  tNlt 
tract  wool#  oft#*  000*90  at  IN#  free  torfac#  of  IN#  avtal  00  cHNer  tide 
of  IN#  torfac#  oolte#  tin#.  So#  cractt  fro*  tN#  toNttrate  were  alt# 
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The  constraints  imposed  by  capillarity  on  the  disposal  of  the  heat 
released  in  crystal  growth  hare  been  evaluated  in  several  recent  publications 
(1-8).  These  evaluations  have  increased  the  accuracy  with  which  the  upper 
Uniting  rates  of  crystal  growth  (i.e. ,  the  heat  transport  limited  rates)  can 
be  computed.  The  rates  of  growth  of  crystals  into  undercooled  metal  melts  are 
thought  to  approach  these  upper  limiting  rates  closely,  reflecting  relatively 
high  rates  of  the  melt  *  crystal  interfacial  rearrangement  processes.  In  this 
near  heat  transport  limited  regime  the  displacement,  aT^,  of  the  crystal-melt 
interface  tcmperaiure,  Tj,  from  the  thermodynamic  equilibrium  temperature,  T  , 
is  negligible  compared  with  aT  «  T^-T,  where  T  is  the  ambient,  i.e.,  heat 
sink,  temperature.  However,  if  the  interface  moves  &T ^  must  be  displaced 
from  zero  by  an  amount  proportional,  near  equilibrium  (Tjn>>aT<)  and  with  the 
fraction  (f)  of  interfacial  growth  sites  constant,  to  the  rate,  V,  of  the 
observed  motion.  In  particular,  (9) 


V  -  V. 


ASf  ATj 


0) 


where  a$^  is  the  molar  entropy  of  fusion  at  T^,  R  is  Avogadro's  number,  and 
kg  is  the  Boltzmann  constant.  VQ  is  a  kinetic  constant  proportional  to  the 
interfacial  rearrangement  frequency  and  f.  It  can  hardly  exceed  the  limit, 
which  should  be  of  the  order  of  the  sound  speed  v  ,  imposed  by  the  collision 
frequency  of  atoms  or  molecules  from  the  liquid  onto  the  interface  (9,  10). 
Thus 


V 


ivs 


(2) 


Exposure  of  the  actual  aT^  and  VQ  should  be  most  favored  under  conditions 
where  the  ambient  undercoolings  and  attendant  crystal  growth  speeds  are  very 
large. 
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Nrhaps  the  highest  measured  crystal  growth  velocities  are  those  of  Ni 
••4  Co  Into  their  pure  undercooled  melts,  reported  by  Walker  (11)  and  Colllgan 
•od  Bayles  (12).  In  their  experiments  columns  of  the  bulk  melt  were  undercooled 
by  varying  amounts  up  to  250*300*  and  then  seeded  at  one  end.  V  was  determined 
from  the  times  at  which  the  freely  moving  dendrite  tip  passed  thermal  sensors 
sited  along  the  column. 

The  growth  rates  at  given  ambient  undercool Ing,  aT,  were  fairly  reproducibl 
to  aT^O.IT^,  but  scattered  widely  at  aT^O.IT^,  presumably  because  of  the 
perturbing  Influence  of  dynamic  nucleatlon  events  which  frequently  accompanied 
the  rapid  growth  in  this  aT  regime.  The  aT  dependence  of  V  was  approximately 
described  by: 

; 

V*A(aT)*  (3' 

with  A  *  0.142  and  0.17  c*/($ec  (K)?)  for  Ni  and  Co,  respectively.  While  still 
on  Its  parabolic  course  V  in  both  metals  had  reached  levels  near  50  meters/sec 
•t  aT  *  O.IT^.  This  velocity  would  require  an  Interfacial  undercooling 

aTj  >  0.0151^  since  vt  In  these  liquids  Is  expected  to  be  -v  3000  meters/sec. 
Therefore  it  should  be  substantially  and  measurably  less  than  the  calculated 
heat  transport  limited  (i.e.,  aTj  *  0)  rate. 

This  paper  presents  the  theory  of  heat  transport  limited  crystal  growth 
as  modified  to  allow  for  the  restraints  imposed  by  the  interface  rearrangement 
process  and  1ts*app1icat1on  to  the  assessment  of  the  relative  importance  of 
heat  transport  and  Interface  kinetics  in  the  growth  of  NI  crystals  in  their 
pure  undcrcoolcd  melts. 


t 
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For  an  Isothermal  dendrite  (paraboloid  of  revolution)  with  tip  radius  R 
propagating  at  constant  velocity  V,  solution  of  the  heat  flow  equation  leads 
to  the  relationship 


[kL/ULLv)]UT  -  AT.)  *  p  exp(p)E-j(p) ,  (4) 

where  p  «  VR/2<L  is  the  Peclet  number,  Ej  is.  the  exponential  integral, 
is  the  thermal  diffusivity  of  the  liquid,  Ly  is  the  latent  heat  per  unit 
volume  and  k^  is  the  liquid  thennal  conductivity.  Equation  (4)  provides  one 
relationship  between  the  three  variables  AT,  V,  and  R.  The  theoretical  work 
of  Langer  and  Muller-Krumbhaar  (1-2)  establishes  that  a  second  relationship 
can  be  obtained  from  morphological  stability  theory.  Further,  an  excellent 
approximation  to  the  more  rigorous  result  of  Langer  and  Muller-Krumbhaar  can 
be  obtained  from  considerations  of  the  morphological  stability  of  a  growing 
sphere.  Such  an  approximate  theory  has  been  shown  to  provide  an  excellent  fit 
to  the  extensive  dendritic  growth  data  on  succinonitrile  (5-6).  We  will  use 

previous  results  on  the  effect  of  interface  kinetics  on  the  morphological 

.  • 

stability  of  a  sphere  to  obtain  a  second  relationship  between  aT,  V,  and  R. 

By  combining  this  second  relationship  with  eq.  (4),  we  obtain  the  dendrite 
growth  velocity  V  as  a  function  of  bath  undercooling  aT.  .* 

We  consider  growth  into  a  undercooled  liquid  of  a  perturbed  sphere  of 
radius  R  .  Above  a  certain  radius,  perturbations  of  the  form  of  spherical 
harmonics  grow  more  rapidly  than  the  sphere.  From  eq.  (13)  of  reference 
13,  this  radius  (above  which  the  sphere  is  unstable)  is  given  by 

i 

» 

*s  ■  (RV2H1  ♦  (M/N)(L  +  a,)}  X  ' 

(1  +  [1  ♦  4  (M/N) (L  -  !)«,/<!  +  (M/N)(L  +  a,)}*]*5),  (5) 
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In  these  equations  M  =  1  +  (1  +  k$/kL)£,  N  =  1-2,  L  =  %(£-!) (1+2 ),  a-j  =  k^/(Lvp^R*) , 
r  =  y/Lu»  y  is  the  solid-liquid  surface  free  energy,  k  is  the  solid  thermal 
conductivity  and  p^  =  (V/aT^)  is  the  linear  interface  kinetic  coefficient.  We 
have  used  the  subscript  s  on  V,  R,  and  aT  to  indicate  that  these  quantities 
pertain  to  the  spherical  geometry.  Eqs.  (5-7)  are  based  on  the  use  of  Laplace's 
equation  for  the  thermal  field  which,  strictly  speaking,  restricts  their 
validity  to  small  values  of  !<l(aT  -  aT^)/[k^Lv];  however  the  time  dependent 
^-stability  analysis  of  Wey  et  al.  (14)  indicates  a  larger  range  of  validity  of 
the  time  independent  analysis.  Since  the  radius  of  the  dendrite  or  sphere  is 
more  than  an  order  of  magnitude  greater  than  the  nucleation  radius  R*,  we  have 
neglected  any  dependence  of  the  bath  undercooling  on  surface  tension. 

Eq.  (5)  gives  the  radius  above  which  the  perturbations  grow  more  rapidly 
than  the  sphere  itself  (see  discussion  of  absolute  and  relative  stability  in 
reference  13);  using  this  value  for  the  radius  seems  more  reasonable  than 
using  that  above  which  the  perturbation  growth  rate  is  positive. 

In  using  the  results  of  morphological  stability  theory  for  a  sphere  to 
obtain  a  relationship  between  AT,  V,  and  R  relevant  to  dendritic  growth,  it  is 
necessary  to  decide  the  manner  in  which  the  sphere  (aT  ,  V  and  R  )  models 
dendritic  growth.  The  basic  hypothesis  is  that  the  dendrite  radius  is  that 
for  which  the  sphere  is  at  the  stability-instability  demarcation  and  we  take  R 
■  Rs>  Since  V  is  an  interface  quantity  while  AT  depends  in  part  on  the  behavior 
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far  from  the  interface,  it  is  perhaps  more  reasonable  to  require  that  V  *  V 
rather  than  AT  =  aT.  We  will  actually  consider  both  of  these  possibilities. 
We  denote  as  the  RV  method  the  case  where  R  =  R„  and  V  =  V„  and  as  the  RT 
method  the  case  where  R  =  R$  and  aT  =  aTs. 

The  only  adjustable  parameter  is  the  integer  Z  which  is  related  to  the 
form  of  the  perturbation  of  the  sphere.  Since  Z  =  1  and  Z  =  2  perturbations 
grow  more  slowly  than  the  sphere,  we  require  Z £  3.  Presumably,  Z  should 
reflect  the  symmetry  of  the  crystal,  which  for  a  cubic  material  suggests  Z  =  4 
or  6  although  Langer  (2-3)  has  argued  for  Z  =  5  for  a  cubic  material.  For 
Z  ~  4  and  kL  =  k$,  eq.  (5)  reduces  to 

I  • 

R$  =  {R*/4}{83  +  9a^  }{1  +  [1  +•  576^/(83  +  9a])2'\li} 

When  interface  kinetics  are  sufficiently  rapid,  i.e.,  -►  0,  Rs  is  proportion 

to  R*  with  the  proportional ity  constant  dependent  on  Z  and  ks/kL>  For  example 
for  a-j  =  0,  the  above  equation  gives  R&  *  83R*/2.  We  also  note  that  when 
interface  kinetics  are  included  the  dendrite  can  not  be  isothermal  and  eq.  (4) 
which  is  based  on  an  isothermal  paraboloid  of  revolution  is  a  reasonable  but 
not  exact  approximation.  Approximate  corrections  to  eq.  (4)  have  been  derived 
(8). 

In  the  RV  method  R  =  R  and  V  =  V  ,  and  we  use  eq.  (7)  to  eliminate  AT 
from  eq.  (5)  and  write  the  right  hand  side  of  eq.  (5)  as  a  function  of  R  and 
V.  For  a  finite  rate  of  the  interface  process  eq.  (5)  gives  a  transcendental 
relationship  between  R  and  V  which  is  solved  numerically.  Eq.  (4)  with 
ATi  *  V/uj  is  then  solved  numerically  to  obtain  V  as  a  function  of  AT.  For 
infinite  interface  kinetics  (i.e.,  v$  =  •»),  eq.  (5)  can  be  solved  explicitly, 


# 


R2  *  [2TMrkL/(LvV)][1  +  ML/N]. 

In  the  RT  method  R  =  R  and  aT  *  AT  ,  and  for  a  given  aT,  R  can  be  calculated 

« 

frotn  eq.  (5);  eq.  (4)  is  then  solved  numerically  to  find  V  as  a  function  of 

AT. 

Calculations  have  been  carried  out  for  the  growth  of  nickel  dendrites. 

The  following  values  of  the  physical  properties  were  used: 


k^  =  k$  =  0.9  J/(cm  s  K) 

=  0.155  cm2/s 
Lv  =  2350  J/cm3 
Y  =  0.4  (I0'i‘)  J/cm2 
'  TMr  =  2.9  (10*5)  K  cm 
v$  »  3.0  (10s)  cm/s 
vi  *  (VsASf/NkBTM)  =  2.0  (102)  cm/sK 

The  values  of  L  ,  Tm  and  the  heat  capacity  of  the  melt  were  taken  from 

the  tabulation  of  Hultgren  et  al.  (15)  while  the  density  of  Ni  at  T  was 

m 

3 

set  equal  to  7.9  gm/cm  (16).  We  found  no  report  of  the  thermal  conductivity 
of  molten  Ni  but  the  listed  k$  of  crystalline  Ni  was  obtained' by  extrapolating 
the  results  in  the  tables  of  Touloukian  et  al.  (17)  to  T  .  Since  these  authors 
presented  arguments  that  the  change  in  the  thermal  conductivity  of  Fe  upon 
melting  is  quite  small  we  have  assumed  that  k^  =  k&  for  Ni.  The  reduced, 
or  scaled,  crystal  melt  interfacial  tension,  o,  is  defined  by 


(  S  )  1/3  r1 

V  uv 
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where?  Is  the  molar  volume.  Studies  of  the  kinetics  of  crystal  nucleation 

in  undercooled  melts  suggest  (9)  that  for  most  metals  a  >  0.6  at  the  onset 

temperature  of  measurable  nucleation  while  the  theory  of  Spaepen  (18)  gives 

an  upper  limiting  value  a  =  0.86  at  Tm  for  metals  crystallizing  to  a  f.c.c. 

-4  2 

structure.  For  our  calculations  we  have  set  y  -  0.4  x  10  J/cm  corresponding 
to  the  assignment  a  =  0.7  to  Ni.  The  speed  of  sound  in  molten  Ni  could  be 
computed  from  [see  Kleppa  (19)]  the  relation 

»,*(»,*>)* 

where  B.  is  the  adiabatic  bulk  modulus  and  p  is  the  density  of  the  melt. 

a 

We  know  of  no  measurement  of  Bfl  of  molten  Ni.  Bundy  and  Strong  (20)  report 

that  the  isothermal  bulk -modulus.  By,  of  crystalline  Ni  at  TR  is  approximately 

1.1  x  10  dynes/cm  and  that  Ni  expands  about  6%  on  melting.  Invoking 

the  Gruneisen  relation  between  B  and  volume  ^  -  yg)  with  Yq  ^  2  we 
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estimate  By  0.85  x  10  dynes/cm  for  molten  Ni  at  TM>  It  is  likely  that  Bfi 
exceeds  By  by  10  to  15%;  thus  our  assignment  v$  'v  3  x  10  cm/s  seems  reasonable. 

The  results  of  the  calculations  are  indicated  by  the  solid  curves  in 
Figures  1-3  for  bath  undercoolings  of  10-200  K.  We  have  indicated  the  ex¬ 
perimental  data  by  the  dashed  line  corresponding  to  Walker's  measurements 
(11),  viz.,  V  =  0.142  (aT)2  in  the  undercooling  range  20-175  K.  While  data  is 

available  up  to  250  K,  the  scatter  becomes  very  large  above  175  K  (11). 

1.8 

Colligan  and  Bayles  (12)  fit  their  data  with  V  =  0.28  (aT)  which  yields 
higher  velocities  for  aT<30K  and  smaller  velocities  for  aT>30K  than  the  quadratic 
expression.  The  maximum  discrepancy  is  a  factor  of  1.4  at  AT  =  175  K  which 
gives  an  indication  of  the  scatter  in  the  measurements.  The  calculated  curves 
In  Figs.  1  and  3  are  based  on  the  RV  method  while  Fig.  2  is  based  on  the  RT 


Method.  Figures'!  and  2  give  results  for  three  different  kinetic  coefficients, 
viz.,  pj  »  »,  200,  and  40  cm/(sK).  The  intermediate  value  of  200  cm/(sK)  was 
estimated  from  the  speed  of  sound  while  the  other  two  values  were  used  to 
indicate  the  sensitivity  of  the  curves  to  the  interface  kinetic  coefficient. 

In  Figs.  1  and  2,  t  »  4  was  used;  Fig.  3  shows  how  the  value  of  t  affects  the 
calculated  results.  The  difference  between  l  -  4,  5,  and  6  is  clearly  smaller 
than  the  discrepancy  between  theory  and  experiment  and  the  scatter  in  the 
experimental  data.  For  m  *  200  cm/(sK),  good  agreement  between  calculated 
and  experimental  values  are  obtained  for  large  values  of  AT.  The  largest 
disagreement  occurs  at  small  values  of  aT»  e.g.,  at  aT  s  20K,  the  experimental 
velocity  (calculated  from  V  -  0.142aT2)  is  57  cm/s  while  the  RV  method  (Fig. 

i 

1)  gives  16  cm/s  and  the  RT  method  gives  30  cm/s.  While  convective  effects, 
which  have  been  ignored'in  calculations,  could  cause  enhanced  growth  rates,  an 
estimate  of  the  type  developed  by  Glicksman  and  Huang  (21)  indicates  that  such 
enhancement  would  be  insignificant  for  the  range  of  supercoolings  considered. 
Better  agreement  betv/een  theory  and  experiment  could  be  obtained  by  using 
different  values  of  the  physical  properties.  For  example,  if  the  thermal 
conductivities  and  diffusivities  are  doubled  or  the  surface  tension  is  halved, 
the  velocities  at  aT  =  20  and  200  K  increase  by  90  and  20%,  respectively. 

It  is  of  interest  to  examine  the  magnitude  of  various  quantities  at  the 
higher  undercoolings.  For  example,  for  the  RV  method  with  Z-  4  and 
AT  *  175  K,  the  calculated  velocity  is  5.14(103)  cm/s  and  the  tip  radius  is 
1.83  (10“5)  cm.  Hence,  the  kinetic  undercooling  AT^  a  25.7  K,  the  dimensionless 
quantity,  kL(AT-AT.)/[KLly]  =  0.37,  and  the  Peclet  number  equals  0.30.  The 
Gibbs-Thomson  lowering  of  the  melting  point  2TMr/R  =  3.2  K,  which  is  negligible 
compared  to  the  total  undercooling  of  175  K.  For  the  RT  method  with  Z  c  4 
and  aT  *=  175  K,  the  calculated  velocity  is  4.28(103)  cm/s  and  the  tip  radius 
Is  2.34(10"5)  cm. 


Our  calculations  show  that  at  the  largest  undercoolings  the  growth  of  Ni 

dendrites  into  their  melt  Is  significantly  retarded  by  the  resistance  of  the 

Interfacial  rearrangement  process.  For  example,  at  an  ambient  undercooling 

«v  0.1  T_  we  calculate  that  the  Interfacial  undercooling  is  -v  0.015  T  or  15S  of 
m  m 

the  ambient.  At  this  large  undercooling  the  measured  speed  of  the  dendrite 
tip  is.  within  the  experimental  and  calculatlonal  uncertainty,  equal  to  that 
calculated  by  assuming  that  the  rate  of  the  interface  process  is  limited  only 
by  the  collision  frequency  of  melt  atoms  with  the  interface.  There  is  no 
evidence  for  thermal  activation  of  this  process. 
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Figure  Captions 

Fig.  1.  Dendritic  growth  velocity  of  nickel  as  a  function  of  bath  undercooling. 

Solid  curves  are  calculated  for  l  »  4  for  three  values  of  the  Interface 
kinetic  coefficient.  Experimental  growth  velocities  are  Indicated  by 
dashed  line  for  which  V  *  0.142  (aT)?. 

Fig.  2.  Dendritic  growth  velocity  of  nickel  as  a  function  of  bath  undercool  Ing 
with  the  same  parameters  as  Fig.  1.  The  assumptions  of  the  RV  method 
and  of  the  RT  method  are  used  in  the  calculations  of  Fig.  1  and  this 
figure,  respectively. 

Fig.  3.  Dendritic  growth  velocity  of  nickel  as  a  function  of  bath  undercooling. 

Solid  curves  are  calculated  for  l  «  4,  5,  and  6  for  the  estimated  kinetic 
coefficient  of  200  ca/(sK).  The  dashed  line  indicates  experimental  data. 
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